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Abstract 



We present a study of three decay modes useful for time-dependent CP 
asymmetry measurements. From a sample of 9.7 x 10^ BB meson pairs 
collected with the CLEO detector, we have reconstructed B^ J/ij^K^, 
XciK^s-- ~^ J/V'7r° decays. The latter two decay modes have 

been observed for the first time. We describe a — >• tt^tt" detection tech- 
nique and its application to the reconstruction of the decay B^ — >■ J/tp Kg. 
Combining the results obtained using vr+vr^ and Kg vr^vr^ decays, we 

determine B{B^ — > J/'0 K^) = (9.5 ± 0.8 ± 0.6) x 10"'*, where the first uncer- 
tainty is statistical and the second one is systematic. We also obtain B{B^ — > 
Xci K^) = (3.9l}|±0.4)xl0-^ and^(S° ^ J/^-k^) = (2.51^;^ ±0.2) x 10"^ 
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CP violation arises naturally in the Standard Model with three quark generations [Q; 
however, it still remains one of the least experimentally constrained sectors of the Standard 
Model. Measurements of time-dependent rate asymmetries in the decays of neutral B mesons 
will provide an important test of the Standard Model mechanism for CP violation 0. 

In this Article, we present a study of J/ipKg, B^ XciKg, and 5° — J/ipir^ 

decays. The latter two decay modes have been observed for the first time. We describe 
a Kg 7r°7r° detection technique and its application to the reconstruction of the decay 

The measurement of the CP asymmetry in B^{^) — > J/ip Kg decays probes the relative 
weak phase between the B^ — 5° mixing amplitude and the b ccs decay amplitude In 
the Standard Model this measurement determines sin 2/5, where /3 = Arg (— V^^V^^/V^^VJ^). 
A measurement of sin 2f3 with B^{B ) xa decays is as theoretically clean as one with 

For the purposes of CP violation measurements, the 5° J/ipir^ decay is similar to 
B^ — >■ D~: both decays are governed by the h ccd quark transition, and both final 
states are CP eigenstates of the same CP sign. A recent search for the 5° D~ 
decay at CLEO established an upper limit on B{B^ D^) |Q. If the penguin (6 dec) 

amplitude is negligible compared to the tree (6 — > ccd) amplitude, then the measurement the 
CP asymmetry in B^{B^) — > J/ipn^ decays allows a theoretically clean extraction of sin 2/5. 
The asymmetries measured with J/ ip K^ and J/ ip 7r° final states should have exactly the 
same absolute values but opposite signs, thus providing a useful check for charge-correlated 
systematic bias in 5-fiavor tagging. If the ratio of penguin to tree amplitudes is not too 
small 0, then comparison of the measured asymmetries in J/ipKg and J/tpn^ modes may 
allow a resolution of one of the two discrete ambiguities (/? — > /3 + vr) remaining after a sin 2(3 
measurement 0. 

The data were collected at the Cornell Electron Storage Ring (CESR) with two config- 
urations of the CLEO detector called CLEO II and CLEO II. V |[. The components of 
the CLEO detector most relevant to this analysis are the charged particle tracking system, 
the Csl electromagnetic calorimeter, and the muon chambers. In CLEO II the momenta 
of charged particles are measured in a tracking system consisting of a 6-layer straw tube 
chamber, a 10-layer precision drift chamber, and a 51-layer main drift chamber, all operat- 
ing inside a 1.5 T solenoidal magnet. The main drift chamber also provides a measurement 
of the specific ionization, dE/dx, used for particle identification. For CLEO II. V, the straw 
tube chamber was replaced with a 3-layer silicon vertex detector, and the gas in the main 
drift chamber was changed from an argon-ethane to a helium-propane mixture. The muon 
chambers consist of proportional counters placed at increasing depth in the steel absorber. 
We use 9.2 fb~^ of e~^e~ data taken at the T{4S) resonance and 4.6 fb"^ taken 60 MeV below 
the T(4S') resonance. Two thirds of the data were collected with the CLEO II. V detector. 
The simulated event samples used in this analysis were generated with a GEANT-based |Q 
simulation of the CLEO detector response and were processed in a similar manner as the 
data. 

We reconstruct both J/ip ^ e~^e~ and J/ip ^ n~ decays and use identical J/ip 
selection criteria for all measurements described in this Article. Electron candidates are 
identified based on the ratio of the track momentum to the associated shower energy in 
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the Csl calorimeter and on the dE/dx measurement. The internal bremsstrahlung in the 
J/ip ^ e~^e~ decay as well as the bremsstrahlung in the detector material produces a long 
radiative tail in the e'^e~ invariant mass distribution and impedes efficient J/ip ^ e^e~ 
detection. We recover some of the bremsstrahlung photons by selecting the photon shower 
with the smallest opening angle with respect to the direction of the track evaluated at the 
interaction point, and then requiring this opening angle to be smaller than 5°. We therefore 
refer to the e~^{'y)e~{'y) invariant mass when we describe the J/ip ^ e^e~ reconstruction. 
For the J /if) ^ fJ'~^fJ'~ reconstruction, one of the muon candidates is required to penetrate the 
steel absorber to a depth greater than 3 nuclear interaction lengths. We relax the absorber 
penetration requirement for the second muon candidate if it is not expected to reach a muon 
chamber either because its energy is too low or because it does not point to a region of the 
detector covered by the muon chambers. For these muon candidates we require the ionization 
signature in the Csl calorimeter to be consistent with that of a muon. 

We extensively use normalized variables, taking advantage of well-understood track and 
photon-shower four-momentum covariance matrices to calculate the expected resolution for 
each combination. The use of normalized variables allows uniform candidate selection criteria 
to be applied to the data collected with the CLEO II and CLEO II. V detector configurations. 
For example, the normalized J/V' —>■ /i"*"/^" mass is defined as [M{fi'^ii~) — Mj/^]/a{M), 
where Mj/^ is the world average value of the J/ip mass ||T0[ and cr(M) is the calculated mass 
resolution for that particular fi~^fi~ combination. The average invariant mass resolution 
is 12 MeV/c^. The normalized mass distributions for the J/ip ^ candidates are shown 
in Fig. |l|. We require the normalized mass to be from —10 to +3 for the J/if) ^ e^e~ and 
from —4 to +3 for the J /tp ^ fJ'~^f^~ candidates. 

Photon candidates for Xci — > J/V' 7 and 77 decays are required to have an energy 

of at least 30 MeV in the barrel region (| cos6'^| < 0.71) and at least 50 MeV in the endcap 
region (0.71 < | cos 6*^1 < 0.95), where 9^ is the angle between the beam axis and the 
candidate photon. To select the vr*^ candidates for ^ J/ip vr*^ reconstruction, we require 
the normalized vr*^ — > 77 mass to be between —5 and +4. The average 77 invariant mass 
resolution for these tt*^ candidates is 7 MeV/c^. We perform a fit constraining the mass of 
each 7r° candidate to the world average value |T^. 

We reconstruct Xci in the Xci ^ J /'>p'~i decay mode. Most of the photons in T(4S') BB 
events come from vr*^ decays. We therefore do not use a photon if it can be paired with another 
photon to produce a vr*^ candidate with the normalized 7r° — > 77 mass between —4 and +3. 
The resolution in the J/ip'^ invariant mass is 8 MeV/c^. We select the xa candidates with 
the normalized xa ^ J/ip'^ mass between —4 and +3 and perform a fit constraining the 
mass of each Xci candidate to the world average value |[To[] . 

The TT+TT^ candidates are selected from pairs of tracks forming well-measured 

displaced vertices. We refit the daughter pion tracks taking into account the position of the 
displaced vertex and constrain them to originate from the measured vertex. The resolution 
in the tt~^tt~ invariant mass is 4 MeV/c^. We select the K^^ 

— > TT+TT candidates with the 

normalized tt+tt^ mass between —4 and +4 and perform a fit constraining the mass 

of each Kg candidate to the world average value ||10 



In order to increase our B^ — > J/ipKg sample, we also reconstruct Kg vr^vr'^ decays. 
The average flight distance for the Kg from B^ — > J/ipKg decay is 9 cm. We find the Kg 
decay vertex using only the calorimeter information and the known position of the e~^e~ 
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interaction point. The flight direction is calculated as the line passing through the e^e~ 
interaction point and the center of energy of the four photon showers in the calorimeter. 
The decay vertex is defined as the point along the flight direction for which the 
product /[M(7i72)] x /[M(7374)] is maximal. In the above expression, M(7i72) and M(7374) 
are the diphoton invariant masses recalculated assuming a particular decay point and 
f{M) is the 7r° mass lineshape obtained from the simulation where we use the known K^g 
decay vertex. For simulated events, the flight distance is found without bias with a 
resolution of 5 cm. The uncertainty in the decay vertex position arising from the 
direction approximation is much smaller than the resolution of the flight distance. We 
select the candidates by requiring the reconstructed Kg decay length to be in the range 
from —10 to +60 cm. After the decay vertex is found, we select the tt^tt" 
candidates by requiring —15 < M(77) — M^o < 10 MeV/c^ for both photon pairs. Then 
we perform a kinematic fit simultaneously constraining M(7i72) and M(7374) to the world 
average value of the 7r° mass [|TU[. The resulting mass resolution is 12 MeV/c^. We select 
the 7r°7r*^ candidates with the normalized vr^Tr*^ mass between —3 and +3 and 

perform a fit constraining the mass of each Kg candidate to the world average value [|r^ . The 
Kg — ^ vr'^vr^ detection efficiency is determined from simulation. The systematic uncertainty 
associated with this determination can be reliably estimated by comparing the Kg — >■ tt^tt" 
and Kg — > 7i~^tt~ yields for inclusive Kg candidates in data and in simulated events. 

The candidates are selected by means of two observables. The first observable is the 
difference between the energy of the candidate and the beam energy, AE = E{B^) — 
-Ebeam- The average AE resolution for each decay mode is listed in Table |I|. We use the 
normalized AE for candidate selection and require \ AE\/a{AE) < 3 for B^ J/ipKg and 

~^ Xci Kg candidates with Kg — >■ 7r~^n~ . To account for a low-side AE tail arising from 
the energy leakage in the calorimeter, we require —5 < AE/a{AE) < 3 for B^ J/ipKg 
with K'^ 7r°7r° and -4 < AE/a{AE) < 3 for 5° J/ipn^ candidates. The second 



observable is the beam-constrained B mass, M{B) = y-Eboam ~P^{B), where p{B) is the 
magnitude of the B^ candidate momentum. The resolution in M{B) is dominated by the 
beam energy spread for all the decay modes under study and varies from 2.7 to 3.0 MeV/c^ 
depending on the mode. We use the normalized M{B) for candidate selection and require 
\M{B) - Mb\/(t{M) < 3, where Mb is the nominal E° meson mass. The AE vs. M{B) 
distributions together with the projections on the M{B) axis are shown in Fig. ^. The 
number of B^ candidates selected in each decay mode is listed in Table |. 

Backgrounds can be divided into two categories. The first category is the background 
from those exclusive B decays that tend to produce a peak in the signal region of the 
M{B) distribution. We identify these exclusive B decays and estimate their contributions 
to background using simulated events with the normalizations determined from the known 
branching fractions or from our data. The second category is the combinatorial background 
from BB and continuum non-BB events. To estimate the combinatorial background, we 
fit the M{B) distribution in the region from 5.1 to 5.3 GeV/c^. As a consistency check, 
we also estimate the combinatorial background using high-statistics samples of simulated 
T(4S') BB and non-BB continuum events together with the data collected below the 
BB production threshold. The total estimated backgrounds are listed in Table H Below we 
describe the background estimation for each decay channel under study. 
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Background for —>■ J/ipK^ with Kg tt+tt^. Only combinatorial background con- 
tributes, with the total background estimated to be 0.3 ± 0.2 events. 

Background for B^ J/ipKg with n^ir^. The combinatorial background is 



estimated to be 0.5 ± 0.2 events. The other background source is S — > J/ip K* |11], with 
K* Ktt^ or K* KgTT with Kg — tt^ti^. The background from these decays is estimated 
to be 0.6 ± 0.2 events. 

Background for B^ Xci Kg. The combinatorial background is estimated to be 0.5 ±0.3 
events. We estimate the background from B decays to the J/ipK^ir final state from the 
samples of simulated events, with the normalizations obtained from the fits to the M{Kti) 
distributions for B^ J/ip KgTr~^ and B^ J/ip K^tt'^ candidates in data. The background 
from B —>■ J/tp KgTT is estimated to be 0.41 ± 0.07 events, and is dominated hj B —>■ J/ip K* 
decays with K* — * Kgir. We find no evidence for B —>■ Xc2 K production and estimate the 
background from B^ — >■ Xc2 Kg to be 0.01 ± 0.01 events. 

Background for B'^ J/il)!:^ . The combinatorial background is estimated to be 0.4lo;3 
events. The AE resolution is good enough to render negligible the background from any of 
the Cabibbo-allowed B J/ipKn^lX) decays, where at least a kaon mass is missing from 
the energy sum. The background from B^ — >■ J/ijjKg decays with Kg tt^tt^ is estimated 
to be 0.38 ± 0.05 events. We estimate the background from B decays to the J/ipim^ final 
state from the samples of simulated events, with the normalizations obtained from the fits 
to the M['ktt) distributions for B^ —>■ J/tpir^Tr^ and B^ —>■ J/ipTT^TT^ candidates in data. 
The background from B J/ipTT^ir is estimated to be 0.2 ± 0.2 events, and is dominated 
by B^ J/ip p'^ decays. 

We use the Feldman-Cousins approach [|T2| to assign the 68% CL. intervals for the 
signal mean for the three low-statistics decay modes {B^ XaKg, B^ — > Jjip-n^^ and 
50 ^ JI^K% with K% ttOttO). We assume ;B(T(4^) ^ E^E^) = i3(T(4^) ^ B^B') 
for all branching fractions in this Article. We use the following branching fractions for 
the secondary decays: B{Jl^ = (5.894 ± 0.086)% [|13l, B{xc\ J li'l) = (27.3 ± 

1.6)% fig, B{Kl TT+TT-) = (68.61±0.28)% fig, andB{K'^g ttV) = (31.39±0.28)% [0. 
The reconstruction efficiencies are determined from simulation. The resulting branching 
fractions are listed in Table |. Combining the results for the two Kg modes used in 5° 
J/ ip Kg reconstruction and taking into account correlated systematic uncertainties, we obtain 
B{B^ J/^K^) = (9.5 ± 0.8 ± 0.6) x 10^^ The measurements of B{B° J/ipK°), 
B{B^ XciK^), and B{B^ — > J/ipn'^) reported in this Article supersede the previous 
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The systematic uncertainties in the branching fraction measurements include contribu- 
tions from the uncertainty in the number of BB pairs (2%), tracking efficiencies (1% per 
charged track), photon detection efficiency (2.5%), lepton detection efficiency (3% per lep- 
ton). Kg — * TT+TT^ finding efficiency (2%), Kg — > vr^vr'' finding efficiency (5%), background 
subtraction (0.01 — 5.5%, see Table D, statistics of the simulated event samples (0.6 — 1.0%), 
and the uncertainties on the branching fractions of secondary decays (see Table |I|). 

In summary, we have studied three B^ decay modes useful for the measurement of sin 2(3. 
We report the first observation and measure branching fractions of the B^ Xci K^ and 
B^ ^ J /ip 7r° decays. We describe a K^ —>■ n^n^ detection technique and its application 
to the reconstruction of the decay B^ J/tpKg. We measure the branching fraction for 
-B" — > J/iIj K^ decays with K^ mesons reconstructed in both tt+tt" and 7r°7r° decay modes. 
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TABLE I. Number of signal candidates, estimated background, average AE resolution, product 
of secondary branching fractions (Bs), detection efficiency, and measured branching fraction. Row 
1 contains the combined value of }3{B^ J/ip K^), rows 2 and 3 contain the individual results for 
the two Kg decay modes. 



Decay 




Signal 


Total 


a{AE) 




Efficiency 


Branching 


mode 




candidates 


background 


(MeV) 


(%) 


(%) 


fraction (xlO~^) 


B^^ 


J/tp K° 












9.5 ±0.8 ±0.6 




TT+TT" 


142 


0.3 ±0.2 


11 


4.04 ± 0.06 


37.0 ±2.3 


9.8 ±0.8 ±0.7 






22 


1.1 ±0.3 


25^ 


1.85 ±0.03 


13.9 ± l.l'' 


8.4t2-i ± 0.7 




Xci 


9 


0.9 ±0.3 


10 


1.10 ±0.07 


19.2 ± 1.3 


3.9tJ| ±0.4 


B^^ 




10 


1.0 ±0.5 


28^ 


11.8 ±0.2 


31.4 ± 2.2'' 


0.25l[J:J^ ± 0.02 



'^The AE distribution has a low-side tail due to the energy leakage in the calorimeter. 
'^Includes the loss of efficiency due to vr^ e'^e~j decays. 
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FIG. 1. Normalized invariant mass of the (a) J/^lJ — > e~^e~ and (b) J/^ ^ /U^A*^ candidates in 
data. The momentum of the J/ip candidates is required to be less than 2 GeV/c, which is slightly 
above the maximal J/^p momentum in i? — > J/ipn decays. The shaded histogram represents 
the luminosity-scaled data taken 60 MeV below the T(4S') showing the level of background from 
non-BB events. 
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FIG. 2. The vs. 

B° ^ J/iIjKI with (c) ^° ^ XciK'^, and (d) 5° ^ J/^/;7r0 candidates. The 

signal candidates, selected using normalized A^' and M{B) variables, are shown by filled cir- 
cles. Below each AE vs. M{B) plot, we show the projection on the M{B) axis with the AE 
requirement applied. The shaded parts of the histograms represent the candidates that pass the 
\M{B) - Mb\/(j{M) < 3 requirement. 
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